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Search for Quark Contact Interactions in Dijet Angular Distributions in pp Collisions 
at = 7 TeV Measured with the ATLAS Detector 

The ATLAS Collaboration 

(Dated: September 28, 2010) 

Dijet angular distributions from the first LHC pp collisions at center-of-mass energy ^/s = 7 TeV 
have been measured with the ATLAS detector. The dataset used for this analysis represents an 
integrated luminosity of 3.1 pb^^. Dijet x distributions and centrality ratios have been measured 
up to dijet masses of 2.8 TeV, and found to be in good agreement with Standard Model predictions. 
Analysis of the x distributions excludes quark contact interactions with a compositeness scale A 
below 3.4 TeV, at 95% confidence level, significantly exceeding previous limits. 



1. Introduction 

At hadron colliders, most events with large transverse 
momentum (pr) transfer occur when a constituent par- 
ton from one of the incoming hadrons scatters from a 
parton in the other. At high pr, these '2 — >■ 2' scattering 
processes are well described within the Standard Model 
by perturbative Quantum Chromodynamics (QCD), the 
quantum field theory of strong interactions. As each 
high-momentum parton emerges from the collision, the 
subsequent parton shower and hadronization create a col- 
limated jet of particles aligned with the direction of the 
original parton. In most of these collisions, two high-p^ 
jets emerge from the interaction. These 'dijet' events are 
particularly useful for measuring quantities associated 
with the initial interaction, such as the polar scattering 
angle in the two-parton center-of-mass (CM) frame, 6*, 
and the dijet invariant mass, rrijj. Precise tests of QCD 
may be carried out by comparing the theoretical predic- 
tions to the experimental distributions. If discrepancies 
between data and QCD are found to be well beyond ex- 
perimental and theoretical uncertainties, this would in- 
dicate that the QCD description needs improvement, or 
that a new process, not included in the Standard Model, 
has appeared. 

This analysis focuses on dijet angular distributions, 
which have been shown by previous experiments [1-4] 
to be sensitive measures for testing the predictions of 
QCD and searching for new processes. Dijet angular dis- 
tributions are well suited to the analysis of early LHC 
data, since they are little affected by the main system- 
atic uncertainties associated with the jet energy scale 
(JES) and the luminosity. QCD calculations predict 
that high-pT dijet production is dominated by t-channcl 
gluon exchange, leading to angular distributions that are 
peaked at |cos0*| close to 1. By contrast, models of new 
processes characteristically predict angular distributions 
that would be more isotropic than those of QCD. 

This Letter reports on the first search with the ATLAS 
detector for quark contact interactions leading to mod- 
ifications of dijet angular distributions in proton-proton 



(pp) collisions at a center-of-mass energy of y/s = 7 TcV 
at the LHC. The data sample represents an integrated 
luminosity of 3.1 pb~^, recorded in periods of stable colli- 
sions, through August 2010. The two distributions under 
study - dijet x distributions, and dijet centrality ratios - 
have been used repeatedly as benchmark measures, and 
will be described in detail below. 

The highest exclusion limits on quark contact inter- 
actions set by any previous experiment [4], for several 
statistical analyses, ranged from 2.8 to 3.1 TeV at 95% 
CL for the compositeness scale A. 

2. Kinematics and Angular Distributions 

The 9* distribution for 2 — > 2 parton scattering is pre- 
dicted by QCD in the partonic CM frame of reference. 
Event by event, the momentum fraction (Bjorken x) of 
one incoming parton differs from that of the other, caus- 
ing the partonic reference frame to be boosted relative 
to the detector frame by an amount which can be deter- 
mined from the dijet kinematics. A natural variable for 
analysis of parton-parton interactions is therefore the ra- 
pidity, y — ^ln( ^]'^^' ), where E is the energy and p^, the 
z-component of momentum, of the given particle. The 
variable y transforms under Lorentz boosts along the z- 
direction as y^y — yB=y — ta,nh~^ {(3b), where (3b is 
the velocity of the boosted frame, and yB is its rapidity 
boost. 

The ATLAS coordinate system is a right-handed 
Cartesian system with the x-axis pointing to the center of 
the LHC ring, the z-axis following the counter-clockwise 
beam direction, and the y-axis going upwards. The polar 
angle 9 is referred to the z-axis, and (f> is the azimuthal 
angle about the z-axis. 

Rapidity differences are boost invariant, so that under 
Lorentz boosts jets retain their shapes in {y,4>) coordi- 
nates. The pseudorapidity, rj = — ln(tan(|)), approaches 
rapidity in the massless limit and can be used as an ap- 
proximation to rapidity. The variables rj and 4> are em- 
ployed in the reconstruction of jets. 

The variable x, used in the first angular distributions 
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considered in this study, is derived from the rapidities of 
the two jets defining the dijet topology {yi and 1/2) ■ For 
a given scattering angle 9*, the corresponding rapidity 
in the CM frame (in the massless particle limit) is y* = 
^ln( il^'^jcogg-.j )■ The variables y* and ys can be found 
from the rapidities of the two jets using y* = ^{yi — 2/2) 
and ys = 5(2/1+2/2)- Then y* may be used to determine 
the partonic CM angle 9*. Additionally, y* is the basis 
for the definition of x- X ~ cxp(|2/i ~ 2/2I) = cxp(2|2/*|). 

The utility of the x variable becomes apparent when 
making comparisons of angular distributions predicted 
for new processes to those of QCD. In QCD, gluon 
(massless, spin-1) exchange diagrams have approximately 
the same angular dependence as Rutherford scattering: 
d7V/dcos6'* oc l/sin''(6'*/2). Evaluation of dN/dx shows 
that this distribution is constant in x- By contrast, the 
angular distributions characteristic of new processes are 
more isotropic, leading to additional dijet events at low 
X- In QCD, subdominant diagrams also cause x distri- 
butions to rise slightly at low x- 

The other important kinematic variable derivable from 
jet observables is the dijet invariant mass, rrijj^ which 
is also the CM energy of the partonic system. In re- 
construction, rrijj is found from the two jet four- vectors: 
nijj = y/{Ei^ + Ei-^Y - {jP^ +p-'^)^, where E and p are 
the energy and momentum of the jets. Both distributions 
used in this Letter arc binned in this variable. 

The second angular distribution considered is the di- 
jet ccntrality ratio, Rc- For this analysis, the detector is 
divided into two pseudorapidity regions: central and non- 
central. Rc is defined as the ratio of the number of events 
in which the two highest px jets both fall into the central 
region to the number of events in which the two highest 
Pt jets both fall into the non-central region. For the cur- 
rent study, the central region is defined as |77i,2| < 0.7, 
and the non-central region as 0.7 < |?7i,2| < 1-3. Since 
new processes are expected to produce more central ac- 
tivity than QCD, their signal would appear as an increase 
in Rc above some mjj threshold, with the increase being 
directly related to the cross section of the new signal. 

Rc distributions are complementary to x distributions 
by being sensitive to different regions of phase space, x 
distributions are fine measures oi 9* and coarse measures 
of rrijj, while the opposite is true for Rc distributions 
as they can be binned more finely in nijj for the given 
amount of data. This gives Rc distributions greater dis- 
crimination in determining mass scales associated with 
hypothetical signals. Ideally, when a signal is present, 
the two distributions together can be used to narrow the 
list of viable hypotheses and to establish the associated 
scale parameters. 

It should be noted that measured x ^^d Rc distri- 
butions are compared to theoretical predictions after jet 
reconstruction in the detector, which smears the parton- 
level distributions. 



3. The ATLAS Detector 

The ATLAS detector [5] covers almost the whole solid 
angle around the collision point with layers of tracking 
detectors, calorimeters, and muon chambers. Jet mea- 
surements depend most strongly on the calorimeter sys- 
tem. The ATLAS calorimeter is segmented in intervals 
of pseudorapidity and (j) to exploit the property that jet 
shapes arc nearly boost invariant in {r],(f>) coordinates. 

Liquid argon (LAr) technology is used in the elec- 
tromagnetic sampling calorimeters, with excellent en- 
ergy and position resolution, to cover the pseudorapid- 
ity range < 3.2. The hadronic calorimctry in the 
range \r]\ < 1.7 is provided by a sampling calorimeter 
made of steel and scintillating tiles. In the end-caps 
(1.5 < \ri\ < 3.2), LAr technology is also used for the 
hadronic calorimeters, matching the outer limits of 
the electromagnetic calorimeters. To complete the 77 cov- 
erage, the LAr forward calorimeters provide both electro- 
magnetic and hadronic energy measurements, extending 
the coverage to \ri\ ~ 4.9. The {ri,(f>) segmentation of the 
calorimeter is sufficiently fine to assure that angular res- 
olution uncertainties in dijet analyses are negligible. In 
the data taking period considered approximately 187,000 
calorimeter cells (98% of the total) were active for event 
reconstruction. 

ATLAS has a three-level trigger system, with the first 
level (LI) being custom built hardware. The two higher 
level triggers (HLT) are realized in software. The HLT 
was not set to reject events accepted by the LI single-jet 
triggers chosen for this analysis. 

4. Event Selection and Reconstruction 

In the current 3.1 pb~^data sample, specific LI jet trig- 
ger selections have been exploited for optimal analysis of 
the angular observables. For both observables, bins of 
nijj are associated with distinct LI jet trigger require- 
ments selected to provide maximal statistics while being 
fully eflacient, as will be detailed for x in Section 7. 

Jets have been reconstructed using the infrared-safe 
anti-fct jet clustering algorithm [6] with the radius pa- 
rameter R — 0.6. The inputs to this algorithm are clus- 
ters of calorimeter cells seeded by energy depositions sig- 
nificantly above the measured noise. Jet four-vectors 
are constructed by the vectorial addition of cell clus- 
ters, treating each cluster as an (E, p) four-vector with 
zero mass. The jet four-vectors are then corrected, as a 
function of and pr, for the effects of hadronic shower 
response and detector material distributions by using a 
calibration scheme based on Monte Carlo (MC) studies 
including full detector simulation, and validated with ex- 
tensive test-beam studies [7] and collision data [8]. 

In order to suppress cosmic-ray and beam-related back- 
grounds, events are required to contain at least one pri- 
mary collision vertex with a position of |z| < 30 cm and 
reconstructed from at least five charged-particle tracks. 
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Events with at least two jets are retained if the highest 
Pt jet (the 'leading' jet) satisfies pi^ > 60 GeV and the 
next-to-leading jet satisfies > 30 GeV. The asym- 
metric thresholds are required so as to avoid suppression 
of events where a third jet has been radiated, while the 
30 GeV threshold ensures that reconstruction is fully ef- 
ficient for both leading jets. Those events containing a 
poorly measured jet [9] with pt > 15 GeV are vetoed to 
avoid cases where such a jet would cause incorrect iden- 
tification of the two leading jets. This criterion results 
in a rejection rate of 0.5% in the current data sample. 
The two leading jets are required to satisfy quality crite- 
ria, such as being associated with in-time energy deposits 
in the calorimeter. To be considered as one of the two 
leading jets, a jet is required to be found within the pseu- 
dorapidity region \t]\ < 2.8, where the jet energy scale is 
known to highest precision. 

5. Monte Carlo Simulations 

The Monte Carlo simulation used for the analy- 
sis presented in this Letter has the following compo- 
nents. The MC samples have been produced with 
the PYTHIA 6.4.21 event generator [10] and the AT- 
LAS MC09 parameter tune [11], using the modified 
leading-order MRST2007 [12] parton distribution func- 
tions (PDF). The generated events are passed through 
the detailed simulation of the ATLAS detector [13], which 
uses GEANT4 [14] for simulation of particle transport, 
interactions, and decays. This yields QCD MC samples 
that have been smeared by detector effects for compari- 
son with collision data. These simulated events are then 
subjected to the same reconstruction process as the data 
to produce dijet angular distributions. 

As the next step, bin-by-bin correction factors (K- 
factors) have been applied to the angular distributions 
derived from MC events to account for next-to-leading 
order (NLO) contributions. The K-factors are defined 
as the ratio N LOme I PYTshow ■ The NLO me sample 
is produced using matrix elements in NLOJET++ [15] 
and the NLO PDF from CTEQ6.6 [16]. The PYTshow 
sample is produced with PYTHIA restricted to leading- 
order (LO) matrix elements and parton showering using 
the modified leading order MRST2007 PDF. Multiply- 
ing the full PYTHIA predictions of angular distributions 
with the bin-wise K-factors, obtained as described above, 
results in a reshaped spectrum which includes corrections 
originating from NLO matrix elements, while retaining 
the non-perturbative effects simulated in PYTHIA. Over 
the full range of X: the K-factors change the normalized 
distributions by up to 6%, with little variability from one 
mass bin to the other. In the case of Rc, the K-factors 
change the distribution by less than 1%. 

The QCD predictions used for comparison to data in 
this Letter arc the end product of the two-step procedure 
described above. 



Other ATLAS jet studies [17] have shown that the use 
of different event generators and different sets of parame- 
ters for non-perturbative effects (e.g. hadronization) has 
a negligible effect on the studied observables in the phase 
space being considered. For the high-pT dijet shape ob- 
servables studied here, x and i?c, differences between 
PDF sets were found to be consistently smaller than the 
uncertainty associated with the CTEQ6.6 PDF set, and 
are not taken into account. 

6. Quark Contact Interaction Term 

The benchmark beyond-the-Standard-Model process 
considered in this Letter is a quark contact interaction, 
which may be used to model the onset of kinematic prop- 
erties that would characterize quark compositcness: the 
hypothesis that quarks are composed of more fundamen- 
tal particles. The model Lagrangian for this benchmark 
process is a four-fermion contact interaction [18-20], the 
analog of the Fermi four-fermion interaction used to de- 
scribe effects of the weak interaction. The effects of the 
contact interaction would be expected to appear below or 
near a characteristic energy scale A. If A is much larger 
than the partonic CM energy, these interactions are sup- 
pressed by inverse powers of A and the quarks would 
appear to be point-like. The dominant effect would then 
come from the lowest dimensional four-fermion interac- 
tions (contact terms). 

While a number of contact terms are possible, the 
Lagrangian in standard use since 1984 [18] is the sin- 
gle (isoscalar) term: Cgggq{A) = |f^^'^7''«'^*^7^*^, 

where /An = 1 and the quark fields '^^ are left-handed. 
The full Lagrangian used for hypothesis testing is then 
the sum of Cqqqq{A) and the QCD Lagrangian. The rel- 
ative phase of these terms is controlled by the interfer- 
ence parameter, ^, which is set for destructive interfer- 
ence (^ ~ +1) in the current analysis. Previous analyses 
[4] showed that the choice of constructive (^ = — 1) or de- 
structive (^ — +1) interference changed exclusion limits 
by - 1%. 

MC samples are calculated in PYTHIA 6.4.21 using 
this Lagrangian, with each sample corresponding to a dis- 
tinct value of A. Angular distributions of these samples 
are processed in the same fashion as QCD distributions, 
including the application of bin-wise K-factors. 

Notably, in addition to quark compositcness, this 
same Lagrangian could be applied to a number of other 
beyond-the-Standard-Model theories (albeit with differ- 
ent coupling constants), so that it serves as a template 
for models of new processes with similar scattering dis- 
tributions. 

7. Kinematic Criteria for Angular Distributions 

The X distributions described here are normalized to 
unit area, (l/A'cv)dA'cv/dx where A'ov is the number of 
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observed events, to reduce the effects of uncertainties as- 
sociated with absolute normahzation. 

Detector resolution effects smear the x distributions, 
causing events to migrate between neighboring bins. This 
effect is reduced by configuring the x bins to match the 
natural segmentation of the calorimeter, by making them 
intervals of constant At;, approximating Ay. This is 
achieved by placing the x bin boundaries at positions 
Xn = e'^" '^^"^ , where n is the index for the lower x bound- 
ary of the n*^ bin, starting with n = 0. In doing this, not 
only is the migration reduced, it is also equalized across 
the span of x- 

The X distributions are divided, in turn, into intervals 
of dijet invariant mass, rrijj. For massless partons, the 
following approximate form shows the dependence of rrijj 
on pT and x- "i^j = ^pri PT2 ■ \/x + 1/x - 2cos(A0). 
Since nijj is the CM energy of the partonic system, new 
processes would be expected to appear in the high mass 
bins. Several rrijj intervals are analyzed to exploit the 
fact that the x distributions in low mass bins would be 
similar to the QCD prediction, while these distributions 
would be modified by new physics processes acting in 
high nijj bins. The sensitivity to these processes depends 
strongly on their cross sections relative to QCD and on 
the number of events in the highest mass bin. 

For X distributions, events are rejected if > 0.75 
or \y*\ > 1.7. The combined criteria limit the rapidity 
range of both jets to \yi,2\ < 2.45. The \y*\ criterion de- 
termines the maximum x, which is 30 for this analysis. 
Taken together, these two criteria define a region within 
the space of accessible yi and y2 where the acceptance is 
almost flat in x- This ensures that the expected shapes 
of the distributions are not significantly changed by the 
acceptance. Also, in low mass bins, the ys criterion em- 
phasizes the contribution from the matrix elements and 
reduces the influence of the effects of PDF convolution. 
In the highest rrijj bin, used for limit setting, the {ysl 
criterion reduces the sample by 16%. These kinematic 
cuts have been optimized through full MC simulation to 
assure high acceptance in all dijet mass bins. 

Since event migration also occurs between bins of rrijj , 
studies of fully simulated jets are used to ensure that mi- 
gration is small. This criterion, along with the require- 
ment of a sufficient number of events, lead to rrijj bin 
boundaries of 340, 520, 800, and 1200 GeV, with no up- 
per bound on the highest bin. As noted earlier, single-jet 
triggers are carefully selected for each bin to be 100% ef- 
ficient. Prescaling of triggers leads to a different effective 
integrated luminosity in each mass bin, with the corre- 
sponding numbers being 0.12, 0.56, 2.0, and 3.1 pb~^ in 
the current data sample for the bins listed above. 

Like the x distributions, the Rc distribution has re- 
duced sensitivity to the absolute JES. However, relative 
differences in jet response in rj could have a significant 
impact on the sensitivity. Hence, for these early studies, 
the T] range is restricted to the more central regions of 



the calorimeter where the JES is uniform to within 1% 
as determined from cross-calibration studies. The Rc re- 
gion has been chosen to end at a maximum of |?7| = 1.3, 
just before the transition region between the central and 
end-cap calorimeters. 

8. Convolution of Systematic Uncertainties 

As mentioned before, the angular distributions have a 
reduced sensitivity to the JES uncertainties compared to 
other dijet measurements. Nevertheless, the JES still 
represents the dominant uncertainty for this analysis. 
The ATLAS JES has been determined by extensive stud- 
ies [21], and its uncertainty has been tabulated in the 
variables ry, pt-, and Ny, the number of vertices in the 
event. The average Ny over the full current data sam- 
ple is 1.7. Typical values of the JES uncertainty in the 
considered phase space are between 5% and 7%. The 
resulting bin-wise uncertainties are up to 9% for the x 
observable, and up to 7% for the Rc observable. 

The dominant sources of theoretical uncertainty are 
NLO QCD renormalization and factorization scales, and 
the PDF uncertainties. The corresponding bin-wise un- 
certainties for normalized x distributions are typically up 
to 3% for the combined NLO QCD scales and 1% for the 
PDF error. 

Convolution of these experimental and theoretical un- 
certainties is done for all angular distributions through 
Monte Carlo pseudo-experiments (PE's). For all events 
in the MC sample 1000 PE's are performed, three ran- 
dom numbers being drawn from a Gaussian distribution 
for each PE. The first is applied to the absolute JES, ob- 
tained from the tabulation described above and assumed 
to be fully correlated across rj. The second number is 
applied to the relative JES, extracted from the same tab- 
ulation, which depends only on 77 and restores the decor- 
relation due to 77 dependence of the energy scale. The 
third number is applied to the PDF uncertainty, pro- 
vided by the CTEQ6.6 PDF error sets. In a fourth and 
final step, the uncertainty due to the NLO QCD renor- 
malization {(J-r) and factorization {fip) scales is found 
by letting fiR and fip vary independently between 0.5, 1 
and 2 times the average transverse momentum of the two 
leading jets, resulting in nine samples drawn from a uni- 
form distribution. In a given PE, the data dijet selection 
criteria described previously are applied. 

9. Comparison of Data to Theory 

In Fig. 1 the measured dijet x distributions are com- 
pared to the QCD predictions, along with 1 a systematic 
error bands determined from the PE's, and statistical 
errors on the data. Figure 2 shows the measured dijet- 
centrality distribution and QCD prediction. The statisti- 
cal uncertainties are obtained using Poisson probabilities. 
In the highest mass bins, the numerator and denominator 
of the ratio typically contain 1 or 2 events each. 
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FIG. 1. The normalized x distributions for 340 < rrijj < 
520 GeV, 520 < nijj < 800 GeV, 800 < mjj < 1200 GeV, and 
rrijj > 1200 GeV, with plotting offsets shown in parentheses. 
Shown are the QCD predictions with systematic uncertainties 
(bands), and data points with statistical uncertainties. The 
prediction for QCD with an added quark contact term with A 
— 3.0 TeV is shown for the highest mass bin rrijj > 1200 GeV. 

To evaluate the agreement between data and QCD in 
Figs. 1 and 2, chi-square goodness-of-fit tests were per- 
formed on each angular distribution under the assump- 
tion that the bin-to-bin correlations are negligible. For 
the X distributions shown in Fig. 1, the chi-square per 
degree of freedom for each dijet mass bin is (from low- 
est to highest) 0.68, 0.83, 0.72, and 0.81, indicating good 
agreement with the QCD prediction. 

Similarly, in Fig. 2 the dijet Rc comparison has a chi- 
square per degree of freedom equal to 0.61, also indicating 
good agreement with the QCD prediction. 

The best fit of the Rc distribution in Fig. 2 is obtained 
for a compositeness scale of 2.9 TeV. This is not statis- 
tically significant, as the QCD prediction lies within the 
shortest 68% confidence interval in 1/A^. 

10. Determination of Exclusion Limits 

Since no signal from new physics processes is appar- 
ent in these distributions, limits have been obtained on 
the compositeness scale A of quark contact interactions, 
based on analyses of the x distributions. The contact 
term hypothesis is tested in the highest dijet mass bin in 
Fig. 1, which begins at rrijj = 1200 GeV. For the x dis- 
tribution in this mass bin, the parameter is defined as 
the ratio of the number of events in the first four x bins 
to the number in all x bins. The upper boundary of the 
fourth bin is at x = 3.32. This choice of the bin bound- 
ary has been determined through a MC study that varies 



FIG. 2. Dijet centraUty ratio, Rc, as a function of rrijj, with 
all events above a mass of 1400 GeV plotted in the last bin. 
Shown are the QCD prediction with systematic uncertainties 
(bands), and data points with statistical uncertainties. The 
prediction for QCD with an added quark contact term with 
A = 2.0 TeV is also shown. 

the number of bins in the numerator, as well as the di- 
jet mass bin, and determines the setting that maximizes 
the sensitivity to quark contact interactions, given the 
current integrated luminosity. 

A frequentist analysis is employed as follows. Predic- 
tions of are obtained for a range of A by interpolation 
between distinct samples generated with different 1/A^ 
values. The QCD sample provides a bound with A — oo, 
and additional samples are generated with A values of 
500, 750, 1000, 1500, and 3000 GeV. A full set of PE's 
is made for each hypothesis to construct one-sided 95% 
confidence level (CL) intervals for F^, and the Neyman 
construction [22] is then applied to obtain a limit on A. 

The result is shown in Fig. 3. The measured value 
of F^ is shown by the dashed horizontal line. The value 
of expected from QCD is the solid horizontal line, 
and the band around it allows one to obtain the 1 a 
variation of the expected limit. The dotted line is the 
95% CL contour of the F^^ prediction for quark contact 
interactions plus QCD, as a function of A and including 
all systematic uncertainties. This contour decreases as a 
function of A since, for a small A scale, there would be 
more events at low x- 

The observed limit on A is 3.4 TeV. This limit is found 
from the point where the 95% CL contour crosses the 
measured F^^ value. All values of A less than this value 
are excluded with 95% confidence. This corresponds to 
a distance scale of ~ 6 • 10^^ fm, from conversion of the 
limit using he. The expected limit, found from the cross- 
ing at the QCD prediction, is 3.5 TeV. 
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FIG. 3. The dashed horizontal Une is the measured F, 
text) and the solid horizontal line is the QCD prediction, with 
a band to illustrate a 1 a variation of the expected limit. The 
dotted curve is the 95% CL exclusion contour for with 
quark contact interactions, used to set the exclusion limit on 
A. 



The impact of systematic uncertainties is as follows. If 
all systematic uncertainties were excluded, the observed 
limit reported above would increase by 6% to 3.6 TeV, 
mainly due to the JES uncertainty. Inclusion of NLO 
scales and PDF uncertainties does not change the limit 
measurably, as shape differences arising from these are 
well below the statistical uncertainties. 

Confirming analyses have been done using a Bayesian 
approach with Poisson likelihoods for x bins, calculated 
using priors flat in 1/A^ or 1/A''. These have resulted in 
observed exclusion limits on A of 3.3 TeV and 3.2 TeV, 
respectively, very close to the limit found in the frequen- 
tist analysis. 

Similarly, an analysis has been performed to establish 
95% CL limits using the dijct centrality ratio shown in 
Fig. 2. The likelihood for Rc is constructed as a prod- 
uct of likelihoods of inner and outer event counts for all 
mass bins, which is then analyzed with a Bayesian ap- 
proach similar to that of the x Bayesian analysis. Using 
priors flat in 1/A^ (1/^^) the observed exclusion limit is 
2.0 TeV (also 2.0 TeV), with an expected limit of 2.6 TeV 
(2.4 TeV), providing an additional benchmark for com- 
parison with other experiments. A weaker limit than the 
one derived from the x analysis is expected due to the 
lower T] acceptance associated with the Rc observable. 

Conclusion 

Dijct angular distributions have been measured by the 
ATLAS experiment over a large angular range and span- 
ning dijct masses up to 2.8 TcV. These distributions arc 



in good agreement with QCD predictions. Using 3.1 pb~^ 
of data, quark contact interactions with a scale A below 
3.4 TeV are excluded at the 95% CL. The sensitivity of 
this analysis extends signiflcantly beyond that of previ- 
ously pubhshed studies. 
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